Giant Monopole Resonance in even-A Cd isotopes, the asymmetry term in nuclear 
incompressibility, and the "softness" of Sn and Cd nuclei 
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The isoscalar giant monopole resonance (ISGMR) in even-A Cd isotopes has been studied by 
inelastic a-scattering at 100 MeV/u and at extremely forward angles, including 0°. The asymmetry 
term in the nuclear incompressibility extracted from the ISGMR in Cd isotopes is found to be 
K T = — 555±75 MeV, confirming the value previously obtained from the Sn isotopes. ISGMR 
strength has been computed in relativistic RPA using NL3 and FSUGold effective interactions. 
Both models significantly overestimate the centroids of the ISGMR strength in the Cd isotopes. 
Combined with other recent theoretical effort, the question of the "softness" of the open-shell nuclei 
in the tin region remains open still. 



The equation of state (EOS) of nuclear matter plays 
an important role in our understanding of a number of 
interesting phenomena such as the collective behavior of 
nucleons in the nuclei, the massive stellar collapse lead- 
ing to a supernova explosion, nuclear properties including 
the neutron-skin thickness of heavy nuclei, and the radii 
of neutron stars [l|, The nuclear incompressibility, 
Kqo, is the curvature of EOS of nuclear matter at sat- 
uration density [3]. is, thus, a measure of nuclear 
stiffness and thereby imposes significant constraints on 
theoretical descriptions of the effective nuclear interac- 
tions. However, even more stringent constraints emerge 
as one studies the evolution of the incompressibility coef- 
ficient as the system becomes neutron rich. Neutron-rich 
systems are sensitive to the poorly-known density depen- 
dence of the symmetry energy and the experiments re- 
ported here are of vital importance in this regard. 

The study of the isoscalar giant monopole resonance 
(ISGMR) provides a direct experimental tool to study 
nuclear incompressibility in finite nuclear systems. The 
centroid energy of ISGMR, Ejsgmr, can be directly re- 
lated to the nuclear incompressibility of finite nuclear 
matter, Ka, as: 



parameterized as p||(| : 



K A 



K vol (l + cA- 1 / 3 ) + K T ((N - Z)/Af + KcouiA-^ 

(2) 

Here, K vo \ is the volume term, directly related to K^, 
c ~ —l [7] t Kcoui is essentially a mo del- independent term 
[1] , and K T is the asymmetry term. Although closely re- 
lated, the finite-nucleus asymmetry term K T should not 
be confused with the corresponding term in infinite nu- 
clear mattcr-a quantity also denoted by K T at times, but 
written here as K^°. Indeed, should never be re- 
garded as the A — y oo limit of the finite-nucleus asymme- 
try K T ;6|. Yet the fact that K T is both experimentally 
accessible and strongly correlated with is vital in 
placing stringent constraints on the density dependence 
of the symmetry energy. Recall that K%° is simply re- 
lated to a few fundamental parameters of the equation of 
state d: 



K™ = K sym - 6L 



L. 



(3) 
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m < r z > 



(1) 



where, m is the nucleon mass and < r 2 > is the mean 
square radius of the nucleus [1, [j| . Ka may be further 
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where Qo the "skewness" parameter of symmetric nuclear 
matter and L and K syul , respectively, are the slope and 
curvature of the symmetry energy. It is the strong sensi- 
tivity of to the density dependence of the symmetry 
energy that makes the present study of critical impor- 
tance in constraining the EOS of neutron-rich matter. 

This asymmetry term, K T , can be studied over a se- 
ries of isotopes for which the neutron-proton asymmetry, 
(N — Z)/A, changes by an appreciable amount. The first 
such investigation was carried out by Sharma et al. f]~oj ] 
and later much improved by Li et al. over the even-^4 
ii2 24gn isotopes [llj, ll2j. Li et al. obtained a value 
K T = — 550±100 MeV, a number consistent with val- 
ues indirectly obtained from several other measurements 



[iiUilHi- Most intriguingly, the ISGMR centroid ener- 
gies of the Sn isotopes were found to be consistently lower 
than those predicted by relativistic and non-relativistic 
calculations by as much as 1 MeV [HI 03. It bears noting 
that the interactions used in these calculations were the 
very same that reproduced the ISGMR centroid energies 
in "standard" nuclei, 208 Pb and 90 Zr, very well, lead- 
ing to the question: Why are the Sn isotopes so fluffy? 

To confirm the value of K T obtained from the Sn iso- 
topes, and to further explore the lack of success of the 
relativistic as well as non-relativistic calculations in cor- 
rectly obtaining the centroids of the ISGMR strength 
in the open-shell nuclei, we have measured the ISGMR 
strength distributions in the even-A !06, 110-116^ - 1SQ _ 
topes. The asymmetry parameter, (N — Z)/A, changes 
by as much as 83% along this chain, making these nu- 
clei very attractive, like the stable Sn isotopes, from the 
point of view of extracting the asymmetry term in the 
nuclear incompressibility. Preliminary results from this 
investigation have been presented previously [13, EH ■ 

The experiment was performed at the Ring Cyclotron 
facility of the Research Center for Nuclear Physics 
(RCNP), Osaka University, Japan. 4 He particles at an 
incident energy of 100 MeV/u were scattered off self- 
supporting even- A 1 6 i 1 io-ii6q ( j targets; highly-enriched 
Cd targets (>93%) with thicknesses ranging from 5 to 6.5 
mg/cm 2 were used. Elastic as well as inelastic scatter- 
ing measurements were performed over a wide range of 
angular settings-elastic scattering over 3.4° — 19° and in- 
elastic scattering at extremely forward angles (0°~9.8°). 
The primary justification for the difficult-to-do extremely 
forward angle measurements lies in the angular distri- 
bution patterns, which exhibit clear distinction between 
various multipoles at these angles and, with the ISGMR 
cross sections peaking at 0°, it is extremely important to 
make a measurement as close to 0° as possible. 

The scattered a particles were momentum analyzed 
by the high-resolution magnetic spectrometer, Grand 
Raiden [19 ]. and focused onto the focal-plane detector 
system comprised of two MWDCs [2(| and two plastic 
scintillator counters. The MWDCs allow measurement 
of the horizontal and vertical coordinates of the impact 
position of the a particle on the focal plane. This, in 
turn, allows the determination of the angle of incidence 
on the focal plane and the momentum of the scattered a 
particles. Using the ray-tracing technique for trajectory 
determination of scattered particles, energy spectra were 
obtained for specific scattering angles by subdividing the 
full angular opening. The focal-plane detector covered 
the excitation-energy range of E x ~ 8 MeV to 31 MeV. 
Energy calibration runs were carried out with a 12 C tar- 
get at every angle for each target. The Grand Raiden 
spectrometer was used in the double-focusing mode in 
order to eliminate the instrumental background [l2l l2l| . 
The background-subtracted spectra at an "average" an- 
gle of 0.7° for various Cd isotopes are shown in Fig. 1. 

Elastic-scattering cross sections were used to extract 
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FIG. 1: Excitation-energy spectra at an "average" angle of 
0.7° for the even-even Cd isotopes investigated in this work. 



the optical model parameters (OMP) for these beam- 
target combinations. The parameters were determined 
using the "hybrid" potential proposed by Satchler and 
Khoa [22j and found to work very well for a-scattering 
at medium energies (see, for example, Refs. [HI, [H, HI ) . 
In this procedure, the density-dependent single-folding 
model with a Gaussian a-nucleon potential was used to 
determine the real part of the optical potential. The 
computer codes SDOLFIN and DOLFIN [H] were used 
to calculate the shape of the real part of the potential and 
the form factors, respectively. For the imaginary term, a 
Woods-Saxon potential was used and its parameters, to- 
gether with the depth of the real part, V, were obtained 
by fitting the elastic-scattering cross sections using the 
minimization of chi-square technique, with the help of 
the computer code PTOLEMY 26]. Using the known 
B(E2) values from the literature [27] and the OMP thus 
obtained, the angular distributions of differential cross 
sections were calculated for the 2+ states. An excellent 
agreement between the calculated and experimental an- 
gular distributions of differential cross sections for the 2f 
states established the appropriateness of the OMP. 

The inelastic-scattering cross sections were sorted into 
1 MeV bins to reduce statistical fluctuations. The 
experimentally-obtained spectra consist of contributions 
from various multipoles. In order to extract the IS- 
GMR contribution from these spectra at different scat- 
tering angles, multipole-decomposition analysis (MDA) 
was performed. The experimental double-differential 
cross sections are expressed as linear combinations of 
calculated distorted-wave Born approximation (DWBA) 
double-differential cross sections for different multipoles 
as follows: 

d 2 a ex n0c. m .,E x ) x ^ , d 2 aS WBA (9 c . m .,E x ) 

d ndE - 2^ aL ^> dndE 

(4) 

where L is the order of the multipole and cil(E x ) is the 
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TABLE I: Lorentzian-fit parameters for the ISGMR strength distributions in the Cd isotopes investigated in this work. Also 
presented are the various moment ratios calculated over the excitation-energy range 10.5 MeV-20.5 MeV, as well as those from 
calculations for the FSUGold, NL3 and SLy5 (without pairing) interactions; mt is the fcth moment of the strength distribution: 
rrik = J S(E X ) dx. For comparison, ISGMR parameters from Lui et al. (Gaussian fits) are also provided, where available 



Target E ISG mr (MeV) T (MeV) Eisgmr (MeV) ^Jm 1 /m- 1 {MeV) y/mz/m^MeV) 









Ref. 


[24] 






Experiment 


FSUGold 


NL3 


SLy5 a 


iUb Cd 


16.50±0.19 


6.14±0.37 






16.06±0.05 


16.83±0.09 


16.27±0.09 


16.73 


17.25 


16.92 


108 Cd 
















16.65 


17.17 




110 Cd 


16.09±0.15 


5.72±0.45 


15.71 


+0.11 
-0.11 


15.72±0.05 


16.53±0.08 


15.94±0.07 


16.59 


17.09 


16.65 


112 Cd 


15.72±0.10 


5.85±0.18 






15.59±0.05 


16.38±0.06 


15.80±0.05 


16.50 


17.00 


16.50 


114 Cd 


15.59±0.20 


6.41±0.64 






15.37±0.08 


16.27±0.09 


15.61±0.08 


16.38 


16.90 


16.47 


116 Cd 


15.43±0.12 


6.51±0.40 


15.17 


+0.12 
-0.11 


15.19±0.06 


16.14±0.07 


15.44±0.06 


16.27 


16.77 


16.36 



mi/mo (MeV) 



a Ref. m 
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FIG. 2: (Color online) The ISGMR strength distributions in 
the Cd isotopes investigated in this work. The solid lines 
represent Lorentzian fits to the data. 



percentage of the energy-weighted sum rule (EWSR) for 
multipolarity L. DWBA cross sections corresponding to 
100% EWSR were calculated using transition densities 
and sum rules provided in Refs. [28|, [29|. DWBA cal- 
culations were performed for up to a maximum angular- 
momentum transfer of L=7; addition of higher angular- 
momentum-transfer terms resulted in minimal to no 
change in the strength distributions. Further details on 
the MDA can be found in Refs. [IJIMSij]. The isovector 
giant dipole resonance (IVGDR) contribution was sub- 
tracted out of the experimental spectra prior to the fit- 
ting procedure. Photonuclear data were used in conjunc- 
tion with DWBA calculations based on the Goldhaber- 
Tellcr model to estimate the IVGDR differential cross 
section as a function of scattering angle [32j . 

The coefficients, ao(E x ), extracted from the MDA were 
used to obtain the ISGMR strength distributions; these 
are shown in Fig. 2, along with Lorentzian fits to the 
data. The choice of Lorentzian shape is arbitrary; the 
final results are not affected in any significant way if, for 
example, a Gaussian fit is used instead. Also, as dis- 



cussed in Refs. [O, [13]) the excess strength at the higher 
excitation energies is attributable to the mimicking of 
the L — angular distribution by components of the nu- 
clear continuum from the three-body channels, such as 
the forward-peaked knock-out process wherein protons 
and neutrons are knocked out by the incoming a projec- 
tiles [ID]. This conjecture is supported by measurements 
of proton decay from the isoscalar dipole resonance (IS- 
GDR) at backward angles: no such spurious strength was 
observed in spectra in coincidence with the decay protons 
[34l HH. The parameters of the Lorentzian fit, and the 
customary moment ratios, y/mi/rri-x, -J m^/mi, and 
mi/mo, computed for consistency from both the exper- 
imental and theoretical strength distributions over the 
energy range 10.5 MeV-20.5 MeV, are provided in Ta- 
ble I. The upper limit of this range has been chosen to 
minimize contributions from the excess strength starting 
at around that energy. It should be clear, however, that 
the comparison between theory and experiment through- 
out this work is for the same range of excitation energy; 
this identical energy range is crucial and decisive in the 
discussion below. 

The moment ratios mi /mo for the Cd isotopes are pre- 
sented in Fig. 3. Also displayed are theoretical results 
extracted from the distributions of isoscalar monopole 
strength computed in a relativistic random phase ap- 
proximation (RP A) using the accurately calibrated NL3 
(^=271 MeV) [1| and FSUGold (^=230 MeV) [H 
effective interactions. A detailed description of the rel- 
ativistic RPA formalism and its implementation may be 
found in Ref. [39|. We note that the FSUGold model em- 
ployed here has been successful in reproducing ISGMR 

centroid energies in 90 Zr, 144 Sm, and 208 Pb [4(|. Yet, this 
same model significantly overestimates the ISGMR ener- 
gies in the Sn isotopes [15[ . Perhaps not surprisingly, our 
present theoretical results overestimate the centroid ener- 
gies in the nearby Cd isotopes as well. The use of the NL3 
effective interaction, with an incompressibility coefficient 
significantly larger than FSUGold, exacerbates the dis- 
crepancy between theory and experiment even further. 
Likewise, in a recently available calculation (36j within 
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FIG. 3: (Color online) Systematics of the moment ra- 
tio, mi/mo for the ISGMR strength distributions in the 
Cd isotopes investigated in this work. The experimen- 
tal results (squares) are compared with relativistic calcula- 
tions performed using the FSUGold (circles) and NL3 in- 
teraction(triangles). Also presented are results from non- 
relativistic calculations performed using the SLy5 parame- 
ter set in the HFB+QRPA formalism with and without the 
mixed pairing interaction (diamonds and stars, respectively) 
[36| . The solid lines are to guide the eye. 



the Skyrme Hartree-Fock+BCS and quasiparticle RPA 
with the SLy5 parameter set (i^oo=230 MeV, which, in- 
cidentally, reproduces the ISGMR in 208 Pb very well), 
the centroids of ISGMR strength distributions in the Cd 
isotopes (also shown in Fig. 3) are, again, significantly 
larger than the experimentally-obtained results. Thus, 
the question originally posed in Refs. [H, HH, H3] 01 "Why 
are the Sn isotopes so Fluffy" extends to the cadmium 
isotopes as well. 

We conclude this theoretical discussion with a brief 
comment on the possible role of superfiuid (or pairing) 
correlations on the softening of the isoscalar monopole re- 
sponse in the Sn and Cd isotopes. This was investigated 
previously by Civitarese et al. 1411 and more recently by 
Li et al. [42[ and Cao et al. [3{j. Indeed, even with 
inclusion of the pairing effects, using a mixed pairing in- 
teraction (3(|, the centroids of the ISGMR remain well 
above the experimental values (see Fig. 3)-the net ef- 
fect appears to be that of lowering the centroid by only 
-100 keV in 106 Cd to a maximum of -240 keV in 116 Cd. 
Thus, the impact of superfiuid correlations on the com- 
pressibility of a fermionic droplet remains an interesting 
open question 1431 to date, and, in spite of significant the- 
oretical effort [36l l42Tj47| , no single approach has been 
able to simultaneously describe the centroid energies in 
90 Zr, 208 Pb, and the Sn/Cd isotopes. The remaining 
challenge, therefore, is not only to describe the distribu- 
tion of monopole strength along the isotopic chain in tin 
and cadmium, but to do so without sacrificing the enor- 
mous success already achieved in reproducing a host of 
ground-state observables and collective modes. 
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FIG. 4: (Color online) The difference K A - K Co uiZ 2 A~ 4/3 in 
the Cd isotopes investigated in this work plotted as a func- 
tion of the asymmetry parameter, (iV — Z)/A. The values 
of Ka have been derived using the customary moment ratio 
*J rni/m-i for the energy of ISGMR, and a value of 5.2±0.7 
MeV has been used for Kcoui @j- The solid line represents a 
quadratic fit to the data. 



As pointed out earlier, the asymmetry term in the nu- 
clear incompressibility can be studied using Eq. 2. It 
should be noted that the first term varies very little over 
the isotopic chains of Sn and Cd for which the neutron- 
proton asymmetry changes by 81% and 83% respectively, 
and that the KqouI is essentially model independent. 
Therefore, an approximately quadratic relation between 
[Ka — KcouiZ 2 A~ i / Z ] and the neutron- proton asymme- 
try, (N — Z)/A, can be applied to fit the experimental 
data. The result of such a quadratic fit is shown in Fig. 4. 
The asymmetry term in the nuclear incompressibility, 
K T , obtained from this fit is — 555±75 MeV; the quoted 
error includes the effect (—20 MeV) of the uncertainty 
in the Kcoui term. This result confirms, and is in ex- 
cellent agreement with, the value K T = — 550±100 MeV 
obtained from the Sn isotopes [III 0~2| ■ This value is also 
consistent with the K T = — 370±120 MeV obtained from 
the analysis of the isotopic transport ratios in medium- 
energy heavy-ion reactions [13J], K T — —5001^" MeV 
obtained from constraints placed by neutron-skin data 
from anti-protonic atoms across the mass table [14j |. and 
K T = — 500±50 MeV obtained from theoretical calcula- 
tions using different Skyrme interactions and relativistic 
mean-field (RMF) Lagrangians Q. 

In summary, we have measured the strength distribu- 
tions of the isoscalar giant monopole resonance in the 
even-A 106 > 110 - 116 Cd isotopes via inelastic scattering of a 
particles at extremely- forward angles, including 0°, with 
the aim to put on further experimental footing the re- 
sults obtained earlier for the ISGMR in the Sn isotopes. 
The centroids of the ISGMR have been calculated in the 
relativistic RPA using the NL3 and FSUGold effective in- 
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teractions. The calculated centroids for the Cd isotopes 
are significantly larger than the experimentally obtained 
values, similar to the results obtained for the Sn isotopes, 
leaving the question of the "softness" of the open-shell 
nuclei open still. The value of — 555±75 MeV for the 
asymmetry term in the nuclear-matter incompressibility, 
K T , extracted from this measurement confirms the value 
obtained from the study of the Sn isotopes. 
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